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^^ Abstract. The epitaxial growth of complex oxide thin films provide three avenues 

•^ to generate unique properties: the ability to influence the 3-dimensional structure of 

c/^ the film, the presence of a surface, and the generation of an interface. In all three 

cases, a clear understanding of the resulting atomic structure is desirable. However, 
determining the full structure of an epitaxial thin film (lattice parameters, space group, 
^ atomic positions, surface reconstructions) on a routine basis is a serious challenge. 

^ In this paper we highlight the remarkable information that can be extracted from 

2 both the Bragg scattering and inelastic multiple scattering events that occur during 

C Reflection High Energy Electron Diffraction. We review some methods to extract 

^^ structural information and show how mature techniques used in other fields can be 

Cl directly applied to the in-situ and real-time diffraction images of a growing film. These 

collection of techniques give access to both the epitaxially influenced 3 dimensional bulk 
structure of the film, and any reconstructions that may happen at the surface. 
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1. Introduction 

The low temperature epitaxial growth of thin films provides a whole new phase space 
for the structural control of materials compared to bulk solid state chemistry processes. 
The ability to affect both kinetic and thermodynamic factors in the growth process with 
low temperatures and epitaxial stabilization provide a rich play ground to manipulate 
and control the structure of materials. There are a wealth of examples in the literature 
of novel epitaxial structures of metallic, semiconducting, and oxide materials. Examples 
include: hep and bcc Cu and Pd^j; GaN with cubic [2j, zinc-blend [3j, and wurtzite[4j 
structures; cubic[5j and hexagonal [6^^ 7] structures of YMnOa; and layered superlattice 
materials, such as SrCu02-BaCu02[8j. 

For correlated oxides, where the electronic and magnetic properties are sensitive to 
small changes in lattice parameter, symmetry, and atom coordination, the ability to grow 
novel structures could have significant technological implications. Yonezawa et a/.[9j 
have grown V2O3 epitaxially on two different substrates and succeeded in changing the 
a lattice parameter by 4%. This reduction in lattice parameter caused the stabilization 
of metallic behavior down to 2K in a system which otherwise would go through a 
metal-insulator transition at 170K. In the case of EuO - a ferromagnetic semiconductor 
currently being explored as a spin injector material for spintronic applications p^ - the 
effects of epitaxy pT]. in conjunction with doping[12j, could provide a significant increase 
to the Curie temperature. 

The perovskite structure - the basic structure of the well studied high temperature 
cuprate superconductors, the colossal magnet oresistive manganites, and the ferroelectric 
titanates - is one case where evidence of subtle, but critical changes in the structure 
due to epitaxy have been clearly documented. For example, SrRuOa, which in its bulk 
form is orthorhombic due to a small distortion of the Ru-0 octahedra compared to the 
simple cubic perovskite structure, becomes monoclinic when grown epitaxially on the 
cubic SrTiOa substrate. p!3l fT4] This subtle change in the structure leads to a uniaxial 
magnetic anisotropy in the epitaxial films instead of the biaxial magnetic anisotropy 
found in the bulk.^j 

The foundation of epitaxial stabilization is the presence of a specific structure on 
a surface to act as a template for the growing film. However, due to the abrupt change 
in coordination of atoms at a surface, the structure of the surface region can differ 
from that of the bulk[15j. As such, the nature - both structural and chemical - of this 
surface region, and our ability to understand how it infiuences an epitaxial interface 
that is formed on it, is of interest. 

The external manipulation of the doping levels of a superconductor via an epitaxial 
interface with a ferroelectric [T6j is one recent example of the potential application of 
epitaxial interfaces. A second is the non-bulk-like charge states at the interface of polar 
films of LaAlOa and LaTiOa on non-polar SrTiOa substrates p!7l [T8] . The properties of 
these polar/non-polar interfaces are reminiscent of the classic semiconductor 2D electron 
gasses, and basic electronic devices are already being fabricated [19]. However, a detailed 
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understanding of the generated interface, and the mechanism responsible for the metaUic 
interface is not yet clear. 

In addition, epitaxial growth also provides a means of purposefully generating highly 
ordered and controlled surfaces. Theoretical work has suggested that the electronic 
and magnetic properties of polar surfaces such as MgO(lll) on Ag(lll), would be 
very different from their bulk materials properties [20j. Due to the varied possibilities 
of reconstruction on a polar surface, such as defect formation or species absorption, 
atomic reconstruction, and even faceting, experimental verification of such work requires 
a detailed view of both the surface chemistry and surface structure. 

In order to better understand the potential of epitaxy and to take full advantage of 
its capacity to generate novel structures, interfaces, and surfaces with unique properties, 
the ability to fully explore the structure (unit cell, space group, atom positions) of 
epitaxial films is a significant part of the puzzle. Ideally structural analysis would be 
available in situ and would provide information not only on the bulk structure of the 
film, but also give access to surface structures. 

Due to the intrinsically highly-ordered nature of epitaxial films, diflFraction is 
well suited to provide this structural information. The direct results of experimental 
diffraction measurements are the unit cell parameters, the space group (or a small 
selection of possibilities) and the intensity data. Full structural knowledge, which would 
include the atomic positions, is not immediately available from standard diffraction 
techniques due to the lack of phase information, although there are many indirect and 
direct methods that can be applied [2 Ij. 

However, using the standard structural determination techniques and tools on thin 
films is a challenging problem. Primarily this is because of the small quantity of the 
material of interest (films are often < 2000 A thick) and the weak interaction of x- 
rays and matter. This is made worse for oxides, due to the low cross-section for low Z 
elements. Furthermore, the two-dimensional nature of the film makes it impossible to 
reach a significant portion of reciprocal space. This shortens the list of possible observed 
and unobserved peaks and their intensities making standard indirect and direct methods 
for atom position determination a serious challenge. 

Synchrotron based x-ray diffraction can help to overcome the weak interaction issue 
by increasing the brilliance of the x-ray beam by a factor of approximately 10^^ over lab- 
based x-ray generation [22j. This increased brilliance allows for the detailed exploration 
of surface-related x-ray diffraction effects via crystal truncation rods (CTR)[23j. CTRs 
are created in reciprocal space, instead of a 3-dimensional reciprocal lattice, due to the 
breaking of the crystal periodicity perpendicular to the surface of the film. The CTRs 
can provide information about the vertical structure of the surface, and can provide a 
clear distinction between bulk structure and that of the surface if there is a change in 
the lateral periodicity of the surface with respect to the bulk. This can occur due to 
surface reconstructions or surface ad atoms, and leads to CTRs that have contributions 
from both the bulk and surface regions, and superstructure rods which are a direct 
consequence of just the surface. Fong and Thompson [24] highlight recent work using 
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synchrotron radiation for in situ studies of bulk and surface properties of perovskite 
ferroelectrics thin films. 

The weak interaction issue with x-ray diffraction studies of thin films can also 
be solved by using electrons as the probe particle. Transmission electron microscopy- 
based diffraction methods can provide a wealth of information on bulk and interface 
structures of thin films (see, for example, [25l |26l [Ej), and can also provide chemical 
information [27j. However, it requires significant sample preparation which is destructive 
to the sample being studied, and does not allow access to surface related effects. 
An in situ technique allowing both the determination of the full crystallographic 
structure (space group, lattice parameters, and atomic positions) and surface related 
crystallographic structures would be a powerful and desirable tool; the ability to obtain 
this type of information in real-time during film growth would be an added bonus. 

The standard tools for in situ structural determination are Low Energy Electron 
Diffraction (LEED) and Reflection High Energy Electron Diffraction (RHEED). LEED 
probes the surface layer of atoms to be studied, and as such it is the primary tool to 
uncover the crystallographic structure (with immediate visualization of the symmetry 
information) of the surface and any surface reconstructions that are present [28j. 
However, the front view geometry of the LEED measurement limits its use as a real-time 
growth monitor. 

The grazing incidence geometry of RHEED makes it well suited to real-time 
monitoring during thin film growth, as the source of electrons and detection scheme 
are physically removed from growth sources and other ancillary growth equipment. 
Minor design modifications even allow its use under poor vacuum conditions, of up 
to 50 Pa[29j. The geometry also provides a high level of surface sensitivity for the 
15-30 keV electrons that are normally used. RHEED equipment is found in the vast 
majority of epitaxial growth chambers, primarily because the basic RHEED equipment, 
its implementation in a growth chamber, and the initial interpretation of the patterns 
are fairly straightforward. 

Due in part to the geometry, and also to the various inelastic scattering effects 
that come along with electron diffraction, there is a wealth of information that can 
be collected and extracted from RHEED based diffraction. This information runs the 
gamut from topological information all the way to full structural information of both the 
bulk and surface region. There are a number of excellent review papers and books that 
cover the many details of RHEED pattern formation, and the wide range of growth and 
topological information that can be accessed from RHEED [30l |3ll [32l [33l [3ll [351136] . 

This paper concentrates on the possible approaches with RHEED to obtain detailed 
structural information. The need for this type of crystallographic information is strongly 
driven by the nature of complex oxide materials - complex in structure and composition 
- where expitaxially driven changes in atomic coordination and small changes in bond 
lengths can have dramatic effects on their properties. In particular, we highlight recent 
techniques that allow rapid data collection so that they can be used during growth 
processes. Furthermore, we are also primarily interested in techniques that can either 
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be interpreted directly from the image, or with fairly straight forward data manipulation 
and simulation tools. 

2. Standard RHEED information 

The geometry of RHEED is shown in figure [T^, along with a typical RHEED image 
in figure ^p. For an ideal fiat surface, the RHEED diffraction image is a gnonomic 
projection of the intersection of the crystal truncation rods (CTRs) and the Ewald 
sphere, and as such is a series of high intensity points on a dark background. These 
diffraction maxima are from single scattering Bragg diffraction events. Mahan[30] 
provides a review of the basic geometrical arguments for RHEED images. On most 
surfaces, however, defects such as step edges or islands generally limit the distance over 
which long range order, and hence coherent interference, are maintained. Broadening 
of the CTRs from these effects leads to streaky RHEED patterns. The form of these 
streaks, and their dependence on azimuthal angle can help distinguish various forms of 
island growth, and rotationally disordered surfaces (see, for example, Ref. [35j). 

The temporal changes in the RHEED images are a rich source of information 
about the growth process. When the diffraction patterns along a high symmetry crystal 
directions are collected during the growth of the film, they can be arranged to highlight 
growth mode changes as a function of film thickness. Figure [2] shows a series of RHEED 
images where the area of the images that contain the first Laue zone have each been 
integrated perpendicular to the shadow edge to generate a line profile, and then plotted 
as function of time. This particular time series is of EuO grown on LaAlOa with an 
epitaxial relationship of EuO(lOO) on LaAlO3(012) with EuO[100]//LaAlO3[l21]. The 
first 18 seconds in figure [2] is prior to the start of growth and shows the bright specular 
spot (00), and the (2l) and (l2) diffraction spots of the the first Laue zone along the [l21] 
direction of LaAlOa. At the zero in time the shutter is opened and growth is commenced. 
This style of RHEED data presentation shows immediately that a disordered interface 
layer grows directly on the substrate, and that as the film thickness increases, a single 
crystal-like growth commences. At times greater than about 20 seconds, the specular 
spot (00), and the (10) and (lO) diffraction spots of the EuO gain intensity. Additionally 
it shows that there is a change in d-spacing perpendicular to the direction of the incoming 
beam of the single crystal-like film compared to the substrate. In this case, the bulk 
d-spacing for EuO of 5.14A is obtained in the film immediately on top of the disordered 
layer, as compared to the 5.35A d-spacing of the LaAlOa substrate. 

RHEED has also been used extensively to track the thickness of materials grown 
in a layer-by-layer fashion. [33j In the simple picture, the order and then disorder of a 
completed and then partially completed growth layer generates an oscillatory behavior in 
the intensity of the specularly refiected beam. These oscillations are seen most clearly at 
low angles of incidence, and often several degrees away from a high symmetry direction 
to minimize the infiuence of multi-scattering events. [37] One of the key applications of 
these thickness oscillations is the calibration of growth rates. They have been used to 
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Figure 1. a) The basic geometrical setup of RHEED with the Ewald sphere 
construction highhghted. b) An exemplary RHEED image with the shadow edge, 
specular spot, Laue zones, and Kikuchi lines identified. 



carefully study the absorption-controlled growth mechanism of oxides, such as EuQ[38j. 
The oscillations also allow a very precise method to control the growth of superlattice 
structures [8j, and to study the diffusion related relaxation of the growth process [39]. 
The growth of the recent SrTiOa/LaAlOa superlattices, with control of layers down to 
half a unit cell is the most recent example of the striking control that is afforded by the 
realtime monitoring of such oscillations [181 SQl E] • 

In addition to the Bragg diffraction that is the primary information in RHEED 
images, there can also be strong patterns generated by the subsequent diffraction of 
inelastically scattered electrons. These patterns, called Kikuchi patterns and discussed 
in detail later, can be quite sharp and contain significant intensity in the case of crystals 
with very high bulk and surface order. The patterns are easily distinguishable from 
the single scattering Bragg diffraction intensity because they move in a continuous and 
coincident manner when the crystal is rotated. As such, they can be used to align the 
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Figure 2. Integrated RHEED-image profile as a function of time of EuO grown on 
LaAlOs, along the [110] direction of the pseudo-cubic unit cell of LaAlOs. At t < 0, 
the profile is that of the LaAlOs substrate. At < t < 20 seconds, which corresponds 
to approximately SA of growth, the RHEED image completely disappears, except for 
some minor intensity in the spectrally reflected spot. At t > 20 seconds a strong 
RHEED pattern reemerges and shows the bulk d-spacing for (100) EuO, which is 
smaller than that of the (110) pseudo-cubic d-spacing of LaAlOs. The right hand side 
of the figure shows profiles from t = — 5 sec and t = 50 sec. 



azimuthal direction of the sample with respect to the incoming electron beam to a high 
accuracy. These Kikuchi patterns are also traditionally used to aid sample alignment in 
Transmission Electron Microscopy. One major complication associated with the Kikuchi 
patterns is that there can be significantly enhanced intensity where the Kikuchi patterns 
and the Bragg diflFraction patterns overlap. This effect is called a surface- wave resonance 
[34] . and can lead to an order of magnitude increase in the intensity. In fact, the RHEED 
patterns from some samples, such as high quality SrTiOa (figure |3]), are so dominated 
by these enhancements that it can be quite hard to clearly see the Bragg diffraction 
pattern. 



3. RHEED as a surface symmetry tool 

Beyond information on surface morphology and time-dependent growth information, 
RHEED also allows the underlying in-plane symmetry of the CTRs to be extracted. 
Collecting RHEED images as the azimuth angle of the sample is varied allows the high 
symmetry crystal directions to be found. Figure |4] shows the RHEED images of a quasi- 
3D lOOOA EuO film grown on LaAlO3(120) at 00, 450, and 26.560 with respect to the 
LaAlOs [l21] direction. These images allow d-spacing to be determined for the (100), 
(110), and (120) planes respectively, and hence the in-plane unit cell lattice parameters of 
the EuO film. Use of the known lattice parameters of the substrate prior to film growth 
allows for the accurate calibration of the RHEED geometry and detection system. 
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Figure 3. RHEED images of SrTiOa. a) Along the (001) zone axis, the 0th zone 
Bragg diffraction is barely visible over the large background intensity, b) 40off the 
(001) zone axis and there is significant enhancement where the Kikuchi pattern and 
the 0th zone single scattering Bragg diffraction overlap. 



The azimuthal angular spread of the high symmetry zone axes aUows the basic in- 
plane symmetry to be determined. Braun et a/. [42] demonstrated one method of more 
directly viewing this symmetry information. By arranging, as a function of azimuth 
angle, a single horizontal slice through the specular spot of each RHEED image, a map 
of intensity versus the azimuthal angle is obtained. This map ideally produces an image 
similar to that of LEED, and as such gives a direct visualization of the symmetry shown 
by the CTRs. Figure [5] shows the a GaAs /3(2 x 4) reconstructed surface from Braunei 
a/, [p. 

When interpreting RHEED with the Ewald sphere construction method, the 
RHEED image is a gnomonic projection of the Ewald sphere onto the flat phosphor 
screen. All high intensity points in the RHEED image indicate locations where the 
Ewald sphere and the CRT's intersect. As such, a simple transformation, relative to the 
shadow edge of the RHEED image, can be applied to each point in the RHEED image 
to map the diflFraction intensity onto the surface of the Ewald sphere in 3-dimensional 
space. By collecting RHEED images as a function of azimuthal angle and transforming 
each image back to the Ewald sphere which is now rotating with respect to the sample's 
reciprocal lattice, a dense intensity map is built up which will reproduce the CTRs 
in 3-dimensions. In scattering language, each point in the RHEED image is a different 
scattering vector (g = kf — ki) with respect to the reciprocal lattice of the sample. As the 
sample is rotated, a large range of scattering vectors are swept out in all 3 dimensions. 
According to the Von Laue condition, intensity will be seen at those scattering vectors 
which match the CTRs. 

When integrating this 3-dimensional intensity map over the direction perpendicular 
to the samples surface, an image can be obtained that displays the crystallographic 
symmetry of the CTRs. This is a much higher signal-to-noise version of the visualization 
presented by Braun et a/.[42j. Some limited structure factor information, such as 
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Figure 4. RHEED images of a quasi-3D EuO film on LaAlOa at (a) 00 to the LaAlOa 
pseudo-cubic (100) direction, (b) 450 , and (c) 26.560. Diffraction from the first Laue 
zone, with expected lateral spacing, can be seen in (c). 



systematic absences, can also be seen. Experimentally, a precession-free rotation of 
the sample provides a stable shadow edge and specular spot, which readily allows this 
transformation to be done in real-time. 

Figure [6^ shows the reciprocal lattice visualization of RHEED data collected over 
an azimuthal angle range of 360° of a Au epitaxial film grown on a mica substrate. 
Au has a cubic unit cell (Fm3m), so the 3-fold symmetry seen in the reciprocal lattice 
visualization immediately shows that the film grew with the close packed (111) plane 
parallel to the mica surface. To further highlight the information that can be quickly 
extracted from these RHEED generated surface reciprocal lattice images, figure[6)3 shows 
the 360° data from a (0001) 0-terminated surface of ZnO {PQ-^mc)^ taken directly after 
heating to 700 °C in vacuum. In this image superstructure rods are clearly seen and 
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Figure 5. Azimuthal scan of a GaAs /3(2 x 4) reconstructed surface. The main 
azimuths and the surface unit ceh are indicated. Reprinted from [42^ 
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Figure 6. A reciprocal lattice visualization from RHEED data of (a) Au deposited on 
a mica substrate, (b) an 0-terminated ZnO (0001) substrate showing a i?30°(V^ x ^/3) 
reconstruction, and (c) a SrTiOa (001) substrate with the reciprocal lattice spacing 
indicated by the grid. 



indicate a i?30°(\/3 x ^/?>) reconstruction is present on this surface. [43] The radially 
directed streaks of intensity which start at the {10/} CTRs are due to the presence of 
Kikuchi lines in the RHEED images. 

We have found that this visualization method only works well if the intensity of 
the Bragg diffraction is significantly larger than that of the Kikuchi patterns. If it is 
not, the intensity from the Bragg diffraction will not be the prominent feature of the 
images. Figure [6b shows the transformed RHEED data of a 180° scan of a SrTiOs 
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(001) substrate. Generally, SrTiOa substrates show very strong Kikuchi patterns, which 
significantly degrade the transformed image. This particular visualization was built from 
RHEED imaged collected on poor quality substrates which minimized the infiuence of 
the Kikuchi patterns. SrTiOa has a cubic (PmSm) perovskite structure with a reciprocal 
lattice spacing of 0.256 A~^, as indicated by the grid on figure |6p. As expected, the 
{01/} CTRs are missing from the image because their structure factor is very small. 

The spectroscopic make-up of the electrons reaching the RHEED screen has been 
explored in some detail by a number of researchers |44l l45l [46] . They have shown that 
the elastic peak is generally well separated from the structured loss distribution, which 
is dominated by surface plasmons, especially at low incident angles[l6]. For the cases of 
Si (111) |46j, SrTiOa (001), GaAs (100), and InP (100)[45] filtering out electrons that 
lost more than 5eV minimizes the infiuence of the inelastically scattered electrons from 
the RHEED image. This energy filtering removes the diffuse background, Kikuchi lines, 
and surface resonance features allowing for a more accurate structure determination 
from the Bragg peaks. Furthermore, superstructure peaks due to possible long range 
ordering of the surface are better detectable and higher quality visualizations of the 
CTRs are generated. 

4. RHEED as a surface structure tool 

RHEED has been a very attractive tool to explore and understand the large range 
of reconstructions present on the surfaces of semiconductors, especially under various 
growth conditions. As discussed in the previous section, the unit cell of the surface 
reconstruction is a fairly straightforward piece of information to obtain from RHEED. 
However, in order to determine the atomic positions of atoms on the surface, information 
must be extracted from the intensity of the RHEED image. 

Conventionally, integrated intensities of various Bragg diffraction spots - often 
chosen so that inelastic scattering effects do not dominate - are collected as the incident 
beam angle is varied in the range of l0 to ^ 50 for a specific azimuthal angle (see, 
for example j47l HHJ [49l [50]). Averages of symmetrically equivalent spots are often 
also included in the data. An alternative approach is to vary the azimuthal angle at 
a specific incident beam angle and integrate the intensity from the specular spot [51j. 
These intensity versus angle plots are then compared with dynamical calculations of the 
intensity of high energy electrons refiected from the surface. The dynamical calculations 
are multi-beam calculations for elastic scattering and are formulated for the case of two- 
dimensional periodicity parallel to a surface [52l [53l [Ml |35j. Calculations of all possible 
atomic configurations consistent with the specific reconstructed surface unit cell need to 
be compared to the experimental data to determine which atomic configuration provides 
the best fit. 

A slightly simplified version of this multi-beam approach is also used, called the 
one-beam method. [551 [35] This method is based on the idea that if an azimuthal angle 
is chosen away from a crystallographic symmetry axis, the incident beam will see each 
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atomic layer as an approximately homogeneous continuum. The variation of the incident 
beam angle will then primarily probe the changes in the crystal potential in the direction 
normal to the surface. Therefore, the specular spot RHEED intensity data can be 
compared to dynamical calculations of models that include only the interlayer distances 
normal to the surface and the atomic densities of those layers as parameters. 

These types of trial-and-error procedures can be very demanding when the model 
parameter space is large, or there is no good initial guess for the surface structure. 
However, recent papers using LEED for surface reconstruction studies have extended 
the possible approaches beyond the conventional dynamical analysis methods by the 
use of the Patterson Function, an application from the kinematical methods of x-ray 
diffraction. [56i [57l l58l [59] The method takes the sum of transforms of LEED I-V spectra 
over multiple incident directions and momenta transfers which minimizes the effects of 
multiple scattering and can produce an artifact-free Patterson function of the surface 
region. 

A Patterson function is a real space function calculated by taking the the Fourier 
transform of measured diffraction intensities in reciprocal space, 

P{^) = ^\Fhki\^exip{-iq-f). 
hki 

It is equivalent to an autocorrelation function of the electron density of a unit cell and 

so provides information about the interatomic vectors between atoms in a unit cell. 

Although the Patterson function approach is not a direct structural determination tool, 

it allows significant insight into the location and arrangement of atoms within the unit 

cell. [23] Following the ideas behind the use of the Patterson function on LEED data - 

large data sets and multiple incident angles - it has recently been applied to RHEED 

data. 

Abukawa et a/. [60] noted that RHEED images taken as a function of azimuthal 
angle and assembled, as previously described, into intensity maps as a function of the 
scattering vector provide access to a large volume of reciprocal space. These scattering 
vectors span the full reciprocal space above the sample, except for a small funnel shaped 
region around the axis of rotation (see figure [7|). This large volume of reciprocal space 
provides intensity data over the CTRs and any SR present, and as such enables the 
possibility of using intensity data to extract atomic positions from the RHEED data. 
Abukawa et a/. [60] presented an initial data set from azimuthal scan RHEED data of a 
Si(lll)-\/3 X ^/3-Ag surface, and used the Patterson function to see changes in atomic 
positions due to a surface structural phase transition at 130K. 

Romanyuk et a/.[62j used the Patterson method to study the formation of iron- 
silicide from two mono-layers of Fe deposited on Si(lll) from azimuthal scan RHEED 
data. To interpret the data they calculated the Patterson function for a set of reference 
structure models of the Fe-Si surface and compared them with the experimental data. 
They were able to discount several models for the surface atomic structure, but they 
found spots in their experimental Patterson function, which might have been artifacts, 
that they could not reproduce with the models. 
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Figure 7. Observable volume in the reciprocal space, where a quadrant has been cut 
to show internal sections, and s, (s^^, Sy^ Sz) is the scattering vector. (Inset) Azimuthal 
rotation of the Ewald sphere. Reprinted from [61^ 



In a second paper, Abukawa et a/.[6l] noted the analogy of the RHEED data 
collected as a function of azimuthal angle with that of the Weissenberg camera method 
for XRD. They also pointed out that using standard kinematic-based x-ray diffraction 
tools, such as the Patterson method, for structural determination with RHEED data 
which clearly includes inelastically scattered structure is problematic. In conventional 
dynamical analysis of RHEED intensities, inelastic scattering issues are minimized by 
carefully choosing the incident beam conditions and which diffraction spots to study. 
With azimuthal scan RHEED data, this is not possible, so instead Abukawa et a/.[61j 
used a high pass filter set to 5.5 eV below the electron beam energy to remove the 
most significant effects of multiple scattering. This enabled a surface crystallographic 
analysis on a Si(lll)-v^x ^/3-Ag surface, via the Patterson function, to produce atomic 
positions in excellent agreement with the known atomic structure. 

The interatomic distances obtained directly via the Patterson function, a kinematic 
diffraction tool, from electron diffraction images provides a very valuable new tool to 
explore the surface atomic structure of epitaxial grown oxides and the surface upon 
which they are grown. Furthermore, the removal of inelastically scattered electrons 
from the RHEED images, by energy filtering, allows quick access to more interpretable 
data over a large area of reciprocal space. In addition, the energy selection of the 
scattered electrons also allows access to basic chemistry information at the surface via 
Auger events. 



5. in situ Kikuchi diffraction 

Removing the inelastic scattering effects from the densely collected RHEED images in 
order to apply kinematic based diffraction methods does, however, remove a significant 
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part of the remaining crystallographic data that can be generated by diffraction in 
the RHEED geometry. In addition to the Bragg reflections of the RHEED beam, 
typical RHEED images of high quahty samples also show Kikuchi lines generated by 
inelastic scattering processes. For many decades Kikuchi patterns have been used 
by electron microscopists to precisely orient samples in 3-dimensional space as they 
contain a representation of the full 3-dimensional crystal structure. Much of the work 
surrounding their use can be directly applied to in situ and real-time pattern collected 
while doing RHEED. This section provides some background on the technique and shows 
its potential to help distinguish subtle crystal structure changes that could be occurring 
due to the epitaxial growth of complex oxides. 

In the 1930's, while studying thin mica fllms using an electron beam, Kikuchi[63] 
observed a background structure, as well as the expected diflFraction peaks (for a 
historical review, see for example, Ref. [64]). This background structure consisted 
of a series of parallel line pairs, later called Kikuchi lines, which were interpreted in 
a "two-event model". The flrst step of the model involved the incoming collimated 
and mono-energetic beam of electrons being diffused in the crystal by an unspecifled 
incoherent scattering process. In essence, this flrst step generates a point source of 
electrons traveling in multiple directions, inside the material. The second step then 
involves standard Bragg diffraction of the diffused electrons from the planes in the 
material. From any point inside the crystal there are two possible angles for which 
Bragg diffraction from a single set of lattice planes can occur. Therefore, each set of 
crystallographic planes generate two diffraction cones with a large opening angle of 
180^ — 26 and an axis perpendicular to the diffracting lattice planes. These cones will 
be seen on a collecting screen as a pair of lines, i.e. the Kikuchi lines (see flgure[8^). 
Zone axes are deflned as the intersection line between a set of crystallographic planes, 
and since the Kikuchi lines are a projection of the crystallographic planes, any point 
where the various Kikuchi line pairs intersect will be a zone axis. 

In reality, the multiple scattering process required to generate the Kikuchi pattern 
leads to more complex contrasts in the image than simple pairs of lines. [65j The 
diffraction pattern is better described as bands with excess or deflcient Kikuchi lines that 
deflne the edge of the band (see flgure[8]3). Deal et a/.[66j studied the energy dependence 
of the Kikuchi patterns and observed that contributions to the Kikuchi patterns arise 
from electrons that have up to about 80% of their incident energy. The Kikuchi band 
contrast is found to be maximum for electrons with 97% of their incoming energy, or 
having lost 450eV of energy with respect to the original 15 KeV beam. 

These Kikuchi patterns (also called Electron Backscatter Diffraction (EBSD) 
patterns) are used extensively in Scanning Electron Microscopes (SEM) and Tunneling 
Electron Microscopes for sample alignment and texture maps based on crystallographic 
information (see, for example, Ref [68j). During RHEED studies, Kikuchi lines are 
often used to demonstrate qualitatively that the long range order of the sample is good. 
They have also been used quantitatively to determine the average crystal potential and 
the crystal misorientation.[33] However, these Kikuchi diffraction patterns contain a 
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Figure 8. a) Scheme of backscatter Kikuchi pattern formation: the primary beam 
enters the sample under the angle a. The primary electrons are incoherently scattered 
with an intensity distribution indicated by the ellipsis, and subsequently coherently 
scattered by the lattice planes p into pairs of Kikuchi cones (1) and (2) with an opening 
angle of 180^ — 29 {0 is the Bragg angle) around the cone axis, b) A typical Kikuchi 
pattern of Si taken with 15 KeV electrons in a field-emission gun Scanning Electron 
Microscope. Reprinted from [67] 



representation of the full three-dimensional structure of the probed material, and can 
therefore be used to determine the space group and potentially the 3-dimensional atomic 
positions of an epitaxially grown film. To make that practical, the Kikuchi patterns need 
to be collected away from the main Bragg diffraction. 

Baba-Kishi[69] first demonstrated the combined collection of RHEED and Kikuchi 
patterns, using two different phosophor screens, within an SEM to study a thin GaAs 
capping layer deposited on GaAlAs. Horio [70] also investigated the combination of 
RHEED and Kikuchi patterns on a Si(lll)-7 x 7 surface. In this case Horio used a 
conventional RHEED setup and added a second hemispherical screen in line with the 
sample normal. This hemispherical screen, which lead to the name astrodome RHEED, 
was used to enable an acceptance angle of 95^ and to minimize the distortion of the 
Kikuchi patterns that occur due to the standard gnomonic projection of the patterns 
when they are collected on a flat screen. Procedures for the the collection and analysis 
of Kikuchi patterns from spherical surfaces have recently been discussed. [71] 

The particular location of Horio's hemispherical screen does not allow for the use 
of this screen during film growth. In figure [9^ we show a slightly altered geometry, 
with a flat screen, that allows the Kikuchi screen to be monitored throughout the film 
growth processes to provide in situ and real-time data collection. The electron beam 
has an angle of incidence on the sample of about 20, as is standard for RHEED. The 
Kikuchi screen is 5 cm away from the sample with an angle of 450 to the sample normal. 
Figure [9b and c show simultaneous RHEED and Kikuchi images collected from a (001) 
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Figure 9. (a) The two phosphor screen geometry to ahow simultaneous RHEED and 
Kikuchi pattern cohection. The simultaneous (b) RHEED pattern and (c) Kikuchi 
pattern of a SrTiOa substrate at 650°C and under a 1 x 10^^ atoms/cm^ sec flux of 
atomic oxygen. 



SrTiOa substrate after heating to 650^0 under atomic oxygen. The strongest Kikuchi 
band contrast is found with the incident angle of the incoming electron beam at around 
150 with this current setup. Zaefferer [67j notes that maximum electron intensity is 
obtained when the angle of incidence of the electron beam and the angle of emission 
of the Kikuchi patters are similar. This suggests that placing a screen to capture the 
Kikuchi patterns with an angle to the sample normal of as close to 900 as possible, 
while still avoiding the main Bragg diffraction peaks, would be advantageous. However, 
as discussed later, the best location for the screen will depend on the particular sets of 
Kikuchi bands that provide access to the most significant structural information. 
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One major concern with using the Kikuchi patterns for the structural analysis of 
thin films is the influence of diflFraction from the substrate on the observed pattern. If 
the fllm is too thin, the electron beam will interact with the substrate as well as with 
the fllm. In general, 30 keV electrons, as used in this work, have a penetration depth of 
roughly 500 A. This is then decreased to about 50 A because of the low angle of incidence 
of the electron beam with the sample. In agreement with this, C.J. Harland et aL[72\ 
determined an information depth of < 100 A, while Kohl[73j calculated an information 
depth of about 50-60 A for the reciprocally-related electron channeling patterns. This 
approximate probe depth of 50-100 A is further experimentally supported by work of 
Baba-Kishi and Dingley[7i] on bulk NiS2. These results imply that a fllm thickness of 
> 200 A will be adequate to avoid the influence of the substrate on the Kikuchi pattern 
formation. 

To gain qualitative information from the Kikuchi patterns requires two important 
pieces of geometrical information: the location of the pattern center {i.e. the electron- 
source point), and the specimen-to-screen distance. Three methods for this are outlined 
in the literature: conic flttingfTSl 176] . the circular mask technique [76]. and the known 
orientation method [77] . Due to the fact that epitaxial growth occurs on a well 
characterized single crystal substrate, the known orientation method is the easiest and 
most accurate method. 

Once the geometrical information is known, the analysis of the Kikuchi pattern 
to obtain space group information [THJ [TU [64] starts with the analysis of individual 
zone axes {i.e. points where at least two sets of Kikuchi lines cross), to deflne their 
respective point group symmetries. The combination of these point group symmetries 
from several zone axes will then allow the determination of the crystallographic point 
group. [64] From this point on, the determination of the space group relies on the ability 
to calculate some rough lattice parameters from interzonal angles from which to simulate 
the Kikuchi diflFraction pattern for comparison with the experimental pattern. These 
procedures have been fully automated in most of the EBSD units attached to SEMs. 

in situ Kikuchi diflFraction of epitaxial thin fllms allows several simpliflcations over 
the general procedure outlined above for analyzing the images. The flrst is that there 
is normally prior knowledge of the expected structure, so that a full ab-inito structure 
determination is not needed. The second is that the diflPr action pattern of the substrate, 
when collected immediately prior to growth, can be directly compared to the as-grown 
fllm to determine changes in structure and to deflne the exact geometry. Finally, the 
azimuthal rotation of the sample to collect RHEED data, also allows the collection of 
the large portion of the Kikuchi pattern that is generated in the hemisphere above the 
sample, thus allowing many more zone axes to be seen. In our experimental geometry 
we can see 2/3 of the full hemisphere, although due to the flat nature of the phosphorus 
screen the edges of our images are geometrically distorted. 



Figure [10] shows a Kikuchi pattern collected after the growth of a SrRuOs fllm 
on SrTiOa (see Ref. [79j and [80] for growth details) while the fllm is still at 650 X^ 
and under a 1 x 10^^ atoms/cm^ sec flux of atomic oxygen. Several zone axes have 
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Figure 10. The Kikuchi pattern of a SrRuOs film grown on a SrTiOs substrate 
at 650°C and under a 1 x 10^^ atoms/cm^ sec flux of atomic oxygen; zone axes are 
indicated. 



been indexed according to the tetragonal unit cell expected at 650 °C.[14J The image 
quality of these particular Kikuchi patterns is of similar quality to many of the original 
patterns collected in SEMs. There has been significant improvements in hardware for 
the collection of these patterns in SEM instruments, and all of those improvements could 
be applied to the in situ collection of images from thin film growth chambers. In fact, 
it is theoretically possible to directly install an EBSD unit, built for an SEM, directly 
into a thin film growth chamber. 

The basic unit cell can be obtained by locating various zone axes on the images, 
and determining their relative geometry, after accounting for any geometrical image 
distortion. Lattice plane spacing, or d-spacing, can then be determined by measuring 
the angular separation of the Kikuchi line pairs, and using the Bragg equation. However, 
as pointed out by Wilkinson [8 Ij, for 20 keV electrons the Bragg angle for low order 
reflections is quite small and the edge of the Kikuchi band is not sharp, which allows 
only a few parts in 100 accuracy. The improvement in this accuracy for higher order 
Kikuchi lines is hindered by their lack of clarity in the image. 

In SEM based Kikuchi pattern measurements, Wilkinson j8l] demonstrated that it is 
possible to obtain sensitivity to strains of the order of 0.02% on thick Sii-a^Ge^^ epitaxial 
layers grown on planar Si substrates. This was done by careful measurements and 
analysis of interzonal angles. This method is limited to strains that distort the shape 
of the unit cell, although that is of primary interest in understanding the tetragonal 
distortion that occurs with epitaxial growth. However, changes in atomic positions may 
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Figure 11. The calculated Kikuchi pattern for tetragonal SrRuOa with the same 



orientation as the experimental image of figure 10 



not be readily detected by these analysis techniques. Therefore we are also interested in 
looking at the differences in the intensity of specific Kikuchi bands as a route to explore 
other epitaxially infiuenced changes in the structure of thin films. 

In order to sort out what to look for in the experimental Kikuchi pattern to help 
distinguish possible changes in the crystal structure of a film due to effects of epitaxy, 
we implemented a simulation based on a pixel by pixel application of the kinematic 
model of diffraction. A kinematic simulation will generate the correct geometry of 
the images, but will not generate correct intensities as Kikuchi pattern formation is a 
multi-scattering phenomena, and therefore requires a full dynamical diffraction model 
to predict intensities. However, Zaefferer[67j argues that the kinematical intensities 
are not too far off the experimental ones, and that with the use of a simple intensity 
correction procedure based on a two-beam diffraction theory, the simple simulation can 
generate qualitatively good results. We have implemented these corrections, as outlined 
in Zaefferer's work, to simulate Kikuchi patterns from various structural models of 



interest. Figure 11 shows a simulation centered on the [133] zone of tetragonal SrRuOa, 
calculated for indices up to 15, which can be compared directly with the experimental 
image in figure [T0| 

An example of a possible epitaxial infiuence on the structure of a grown film is the 
potential buckling of the Cu-0 planes in the high-pressure phase of SrCu02 when grown 
with compressive strain. SrCu02 is the end member (n = oc) of the homologous high- 
pressure series Sr^_iCu^+i02n,[82]. Due to the simplicity of the crystal structure and 
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Figure 12. The calculated Kikuchi patterns, and structures, of the high pressure 
phase of SrCu02 without buckled Cu-0 planes (a), and with buckled Cu-0 planes (b). 
The light and dark lines at the top and bottom of the image are guides to the eye 
of the respective excess and deficient lines that define the edges of the Kikuchi bands 
that change between the two images. 



the planar nature of its Cu and O ordering, it is often called the parent compound of the 
high temperature superconductors [83] and is considered to be a 2-dimensional quantum 
magnet. [84] An extremely large change in the c-axis of this material was seen when it was 
grown by epitaxy with reduced a and h lattice parameters. [85] One possible explanation 
for this anomalous expansion in the c-axis is the buckling of the Cu-0 planes, suggested 
by Vailionis et a/..[86j The proposed buckling of the Cu-0 planes occurs by the motion 
of just the oxygen atoms, and as such it is extremely hard to see in x-ray diffraction. 
However, Kikuchi diffraction is an ideal tool to look at the potential of an epitaxially 
driven change in the dimensionality of the Cu-0 planes. 

Figure [12] shows several changes that are expected in the Kikuchi pattern, centered 
around the (130) zone axis, if the Cu-0 planes were to buckle by the oxygen atoms 
moving in and out of the a-h plane by 10% of the c lattice constant (Figure fl2b) in 



contrast to a planar Cu-0 structure (Figure 12^). In this part of the Kikuchi pattern 
there is an increase in intensity of the scattering from the {224} planes - indexed 
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according to the enlarged unit cell required to describe buckling - due to buckling 
(highlighted in figure [T2|3). Furthermore, there is also significant reduction in intensity 
expected for the (001) and (003) planes, as well as the cubic-indexed {204} planes for 
the buckled structure (highlighted in figure [12^) . 

There are a number of other expected changes in the Kikuchi pattern for the 
buckling of the Cu-0 plane in SrCu02, but primarily in the high-index planes. Although 
high-index planes are present in the image, they not only have inherently a lower 
intensity but are overlaid with the more intense low-index planes, and locations in the 
pattern are not always present that allow them to be seen. However, these simulations 
suggest that even subtle changes in the structure of epitaxially grown films are quite 
distinguishable in RHEED generated Kikuchi patterns. In addition, the ability to collect 
a Kikuchi pattern from the well characterized substrate just prior to the growth of a film 
provides a remarkable opportunity to minimize unknown experimental and geometrical 
effects in the patterns and clarify the significant aspects of the simulations. This affords 
a tremendous advantage to the final analysis and interpretation of the Kikuchi patterns 
collected from the epitaxial film. 

Recent work on multi-beam dynamic models, by Wilkelmann et a/. [87], which 
include improvements in the handling of inelastic scattering processes and the 
description of the excitation process show very good agreement between simulated 
and experimental patterns of GaNjOOOl}. The simulations have reproduced the fine 
structure around zone- axes and higher-order Laue zone rings seen in the experiments. 
These results, along with continual improvements in hardware [88j and software that are 
being applied to EBSD in SEMs, allow a much more accurate and fine-grained approach 
to inspecting and understanding the experimental Kikuchi patterns collected from thin 
epitaxial film to look for subtle changes in the atomic structure. 

6. Conclusion 

The use of high energy electrons at grazing incidence during the growth of epitaxial 
films generates scattered electrons which can be analyzed to obtain an abundance of 
information about the morphology and the crystallographic structure of the film. The 
collection of the standard RHEED patterns at all azimuthal angles of rotation, and 
ideally with energy filtering to isolate just the elastic scattering, can be reconstructed in 
real-time to provide immediate visualization of the symmetry of the crystal truncation 
rods, highlighting any surface reconstructions that may be occurring, as well as the 
basic in-plane lattice parameters. The intensity from this type of data set can also 
be treated with the Patterson function approach of x-ray diffraction to help reveal the 
atomic positions of atoms on the surface. Furthermore, the electrons that undergo 
inelastic scattering to generate Kikuchi patterns, also allow access to the space group 
and atomic positions of the thin epitaxial film. These Kikuchi patterns can also be 
collected in real-time, and simultaneously with the RHEED data on a separate screen. 
These various methods to collect and analyze the scattered electrons from RHEED 
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allow a remarkably complete set of data which help elucidate the full crystallographic 
structure of the potentially novel structures, interfaces, and surfaces that can be formed 
by the influence of epitaxy on the growth of thin fllms. 
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